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Abstract 


Severe Acute and Middle East Respiratory syndrome coronaviruses (SARS-CoV and MERS- 
CoV) encode multifunctional papain-like proteases (PLPs) that have the ability to process the 
viral polyprotein to facilitate RNA replication as well as antagonize the host innate-immune 
response. The latter function involves reversing post-translational modification of cellular 
proteins conjugated with either ubiquitin (Ub) or Ub-like interferon stimulated gene product 15 
(ISG15). Ubiquitin is known to be highly conserved. among eukaryotes but surprisingly ISG15 is 
highly divergent among animals. The ramifications of this sequence divergence to recognition of 
ISG15 by coronaviral papain-like protease at the structural and biochemical levels are poorly 
understood. Therefore, the activity of PLPs from SARS-CoV, MERS-CoV and mouse hepatitis 
virus (MHV) was evaluated against seven ISG15s originating from an assortment of animal 
species susceptible, and not, to certain coronavirus infections. Excitingly, our kinetic, 
thermodynamic and structural analysis revealed an array of different preferences among PLPs. 
Included in these studies is the first insight into a coronoavirus PLP’s interface with ISG15 via 
SARS-CoV PLP in complex with the principle binding domain of human and mouse ISG15s. 
The first X-ray structure of the full-length mouse ISGI5 protein is also reported and highlights a 
unique, twisted-hinge region of ISG15 that is not conserved in human ISGI5 suggesting a 
potential role in differential recognition. Taken together, this new information provides a 
structural and biochemical understanding of the distinct specificities amongst coronavirus PLPs 
observed and addresses a critical gap of how PLPs can interact with ISG15s from a wide variety 


of species. 


Ub, Ubiquitin; ISG15, Interferon stimulated gene product 15; PLpro, Papain-like proteases; PLP, Papain-like proteases (for viruses that have 
more than one PLpro); MHV, Mouse hepatitis virus; CoV, Coronavirus; vOTU, Viral ovarian tumor domain proteases; IFN, Type-I interferon; 
Ubl, Ubiquitin-like domain; NS and Nsps, Nonstructural proteins; USP, Ubiquitin-specific proteases; ITC, isothermal titration calorimetry 
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INTRODUCTION 


Coronaviruses (CoVs) are enveloped, positive-stranded RNA viruses that cause mild to 
severe infections in a wide range of mammals and birds. Specifically, severe acute and Middle 
East respiratory syndrome coronaviruses (SARS-CoV; MERS-CoV) are well-recognized viral 
pathogens that have emerged from different animal reservoirs to cause deadly disease in humans. 
SARS-CoV first emerged in 2002 with a case fatality rate of 10%, claiming the lives of over 800 
people and infecting more than 8000 " *. Ten years later, MERS-CoV emerged with a shocking 
case fatality rate near 35% and has spread to 27 different countries to date 3. The continuing 
threat of MERS-CoV was recently underscored by one of its most recent outbreaks in Republic 
of Korea. This outbreak quickly led to 36 deaths, which has brought the total MERS-CoV global 


deaths to over 600 “>. 


Similar to other positive stranded RNA viruses, CoVs encode two types of cysteine 
proteases, including the papain-like protease (PLP) and 3C-like protease, also known as the main 
protease. Together, these enzymes cleave the viral polyprotein into 16 different nonstructural 
proteins (Nsps 1-16) in order to generate the membrane-bound replicase complex for RNA 
replication 6 7:89 CoVs encode either a single PLP, termed PLpro, or two PLPs that process a 
total of three cleavage sites within the polyprotein ” '© For instance, SARS-CoV and MERS- 
CoV encode a single PLpro, while other CoVs such as mouse hepatitis virus (MHV) encode for 


both the papain-like protease 1 (PLP1) and papain-like protease 2 (PLP2) (Figure 1a). 


Beyond cleaving the viral polyprotein, PLPs have additional activities that promote virus 
replication. The X-ray structure of the first CoV PLP determined from SARS revealed that these 


enzymes resemble the structure of human ubiquitin-specific proteases (USPs) and are thereby 


known as viral USPs, often acting as deubiquitinating enzymes with the ability to remove the 
post-translational modification ubiquitin (Ub) from target proteins *. Some PLPs are also 
deISGylating enzymes with the ability to reverse the post-translational modification of the Ub- 
like protein interferon stimulated gene product 15 (ISG15) from cellular proteins ''. Such 
activities were implicated in SARS-CoV’s suppression of the innate immune responses, 
particularly antagonizing type-I interferon (IFN) signaling and chemokine and cytokine 
production '2: 13 Simultaneous disruption of both activities, either by mutation in MERS CoV 
PLpro’s Ub/Ub-like binding region ', or a destabilizing mutation in MHV PLP’s ubiquitin-like 
domain (UbI) '° has been observed to prevent antagonization of IFN signaling, chemokine and 
cytokine production, as well as viral pathogenesis respectively. Overall, these studies suggest 
that the deubiquitinating and deISGylating activities of PLPs, sometimes packaged together as 


deubiquitinating enzyme (DUB) activity, are a likely contributor to viral pathogenesis. 


Intriguingly, the precise role of each individual activity in promoting PLPs ability to act 
as an IFN antagonist has yet to be precisely defined. So far, several X-ray structures of PLpro 
and PLP2 bound with Ub molecules have been reported ie However, no structure of any 
CoV PLP in complex with an ISG15 molecule has yet been reported. Due to the lack of 
structural detail about the interface between CoV PLPs and ISG15, understanding how CoVs 
PLPs specifically engage with ISGI5s versus Ub has been especially difficult to understand. The 
lack of structural information has led to problems in defining the differences between CoV PLPs’ 
deubiquitinating and deISGylating activities among different CoVs. A further complication in 
our understanding stems from the fact that the sequence identity of ISG15 among mammals is 
low, ranging from 58% to as low as 35% when a broader range of animals is compared (Figure 


1b). In contrast, Ub is highly conserved among eukaryotes. 


The sequence diversity of ISG15 among animals and the impact of this diversity on the 
recognition of these ISG15s by viral USPs has not been well studied despite the potential 
implications. For example, certain CoVs, such as SARS-CoV and MERS-CoV, are known to 
replicate and survive in a broad range of animals including bats, camels, mice, civets, shrews, 


; 10; 17; 18 
badgers, pigs and humans ~ ” 


. In contrast, other CoVs such as MHV only replicate in mice, 
suggesting that some CoV PLPs may have evolved strict specificity for their single host’s 
ISG15'°. The potential for species-species variances in ISG15 have already been shown to be a 
key factor in other viruses including influenza B where it was shown that NS1 is unable to bind 
mouse ISGI5 unlike its non-human primate and human counterparts limiting influenza B virus 


; é Saks - 20; 21; 22 
infection in mice ~’~’ ™. 


Also, biodiversity that occurs between species within ISG15s was 
recently shown to impact the ability of nairovirus viral ovarian tumor domain proteases (vVOTUs) 
to effectively process certain ISG15s **, Overall, many CoV PLPs have been observed to show 
robust deubiquitinating activities; ‘4 however, recent in vivo studies have started to unveil the 


25; 26 


importance of PLPs deISGylating activity . This has led to a need for a better and more 


detailed understanding of the interactions between CoV PLPs and ISG15 at the molecular level. 


Adding to the ambiguity of how CoV PLPs or other deISGylating enzymes interact with 
ISG15, especially when it comes to understanding species-to-species variations, is the lack of 
available full-length ISG15 structures that have been determined. This is especially important 
when the intramolecular arrangement of the two domains of ISG15 is considered. For example, 
the hinge region of ISGI5 shows a significant degree of sequence diversity between the different 
species (Figure 1b). Although several structures of human ISG15 (hISG15) have been reported, 


no complete structure of ISG15 from another species has been resolved. Moreover, the recent 


elucidation of the structure of the C-terminal domain portion of mouse ISG15 (mISG15) 


prompted questions on the impact of ISG15 biodiversity on the overall ISG15 structural fold = 


To address these critical gaps in our understanding of ISG15 recognition by CoV PLPs, 
PLpro from SARS-CoV and MERS-CoV, as well as PLP2 from MHV, were examined for 
differences in their selectivity among ISG15s from various animals. The X-ray crystal structures 
of SARS-CoV PLpro bound to the C-terminal domain of ISGI5 originating from mouse and 
human were determined and analyzed in conjunction with enzyme kinetic and thermodynamic 
data derived from ITC. In addition, the first X-ray structure of a complete non-human ISG15 
structure, mouse ISG15 (mISGI5), was also determined to elucidate potential sequence and 
structural differences that may account for species specificity of CoV PLPs. Together, these 
studies provide significant and new insight into the CoV PLP’s ability to accommodate the 


structural differences not only between Ub and hISGIS5, but also different species ISG15s. 
RESULTS 
Species specific cleavage of proISG15s by CoV PLPs. 


Currently, there is a paucity of available biochemical data on the impact of species-to- 
species sequence variations within ISG15 and the effect that these differences may have on the 
ability of CoV PLPs to recognize and cleave ISG15. Therefore, we employed a recently 
developed assay for deISGylating vOTUs that takes advantage of the ability of a protease or 
DUB to cleave immature ISG15 **. The PLpro enzymes from SARS-CoV and MERS-CoV, as 
well as the enzyme PLP2 from MHV, were purified and assessed for their ability to cleave 
proISGI5s derived from 7 different species including human, sheep, northern tree shrew, 


jackknife fish, mouse, dromedary camel, and vesper bat (Figure 2). Each of these proISG15 


proteins is appended with the pro[SG15 extension from Homo sapiens (Figure 2). ISG15 from 


jackknife fish was included in the analysis to add a more distantly related ISG15 homologue. 


SARS-CoV and MERS-CoV PLpros are both capable of fully processing proISG15 from 
human, mouse, camel and bat within 60 minutes. SARS PLpro is also able to fully process 
proISG15 from sheep and shrew within that same time but it has little to no activity against 
jackfish proISG15. In contrast, MERS PLpro is fully capable of processing jackfish proISG15 
but it processes proISG15 from shrew and sheep poorly. MHV PLP2 on the other hand displays 
a much narrower specificity for ISG15s from different species. It shows little to no cleavage of 
proISG15s derived from human, sheep, camel, or bat sources but does show modest activity for 
the northern tree shrew. MHV PLP2 shows strong activity towards jackfish and mouse 


proISGI1S5, the latter activity being consistent with MHV’s pathogenicity as a murine CoV. 


Species specific affinity of ISG15 for SARS-CoV and MERS-CoV PLpros 


Analysis of the cleavage patterns in Figure 2 suggests that SARS and MERS PLpros are 
more promiscuous in recognizing and cleaving various proISG15s compared to MHV PLP2, 
which may relate to the fact that these human pathogens are capable of replicating in different 
hosts, e.g. bats, camels and the shrew, in contrast to MHV which can only replicate in the mouse. 
Differences in cleaving ability of pro[SGIl5s by SARS and MERS PLpro may be due to 
sequence and structural differences, which may affect binding affinities for the ISG15s. This 
possibility was explored using Isothermal Titration Calorimetry (ITC) to measure the binding 
affinities of different ISG15s missing the P’-sites (Figure 3). As deISGylases natural substrates 
is an isopetide bond, removal of the P’-sites also allows the avoidance of potential confounding 


factors related to CoV PLPs other prominent function of processing the viral polyprotein to 


promote viral replication. This part of the process requires that PLPs recognize and cleave the 
peptide bond after the LXGG sequences within the polyprotein, and the ability of PLPs to cleave 
this bond has been shown to be reliant on residues flanking the peptide cleavage site &? This 
requisite function is not associated with other viral deISGylating proteases such as nairovirus 
vOTUs”’. Therefore, ITC measurements were performed using mature ISG15. Although the 
mature ISG15s resemble the product by the absence of an isopetide bond, binding affinity could 
be assessed independent of any P’-sites contribution associated with using the ISG15-pro-form 
protein. This would also provide a more quantitative understanding of the thermodynamic 
parameters related solely to the species and viral related variances of the protease and ISGI5s 


respectively. 


Intriguingly, both SARS-CoV and MERS-CoV PLpro have stronger affinity toward 
hISGI5 than shISG15 and mISG15. The dissociation constants (Kg) of SARS-CoV and MERS- 
CoV PLpro for hISG15 are 20.5 + 4.5 uM and 59.3 + 4.5 uM, respectively. Unlike their 
affinities for hISG15, the affinities of both PLpros for shISG15 and mISGI5 were significantly 
weaker and as a result, a competitive ITC binding assay (Table 2; Figure 3) had to be used to 
determine their Kq values. In regards to SARS-CoV PLpro, ITC analyses revealed similar Kq’s 
of 200 + 41 uM and 198 + 64 uM for shISG15 and mISG15. More divergence in affinity for 
mISGI15 and shISGI5 were seen with MERS-CoV PLpro. Like SARS-CoV PLpro, affinity of 
MERS-CoV PLpro for shISG15 was similar with a Kg of 147 + 36 uM. However, a substantially 


weaker affinity of MERS-CoV PLpro for mISG15 was observed (Kg of 376 + 53 uM). 


The measured thermodynamic parameters revealed that both SARS-CoV and MERS- 
CoV PLpro follow a similar trend and prefer hISG15 over shISG15 and mISG15. Specifically, 


the association of shISG15 and mISGI5 to PLpros is unfavorable with increasingly higher 


entropic factors. Although PLpros originating from MERS-CoV and SARS-CoV exhibited a 
similar preference for hISG15, the thermodynamics driving the affinity differ. The SARS-CoV 
and hISGI5 binding event was driven by enthalpic factors and was slightly entropically 
unfavorable. In contrast, the MERS-CoV affinity of hISG15 was more balanced possessing an 


entropic and enthalpic component. 


Although some DUB proteases that possess deISGylating activity have been observed to 
engage only the C-terminal domain of ISG15, SARS-CoV has been suggested to interact with 


. 5 15; 23; 27 
both domains! !: 73: 


. To gain insight into the relative contributions of the C-terminal domain 
to their full-length counterparts, ITC was performed on utilizing both mouse and human C- 
terminal variants of ISG15. The Kg of 57.6 $3.21 uM was observed for ChSIG15. Although the 
binding event still was enthalpically favorable, this contribution was four-fold less that that 
observed for hISG15. Conversely, it possessed a markedly favorable entropic component over 
the slightly disfavorable one of its full-length counterpart. Not unsurprisingly, the trend of 
SARS-CoV having a stronger affinity toward hISG15 over mISGI5 was also observed in the C- 
terminal domain of each substrate. The Kg for CmISG15 was considerably weak at 1,870 + 173 
uM showing marked decreases in both entropic and enthalphic binding contributions suggesting 


that similar to hISG15, mISG15’s N-terminal domain may also be involved in the protease- 


ISG15 binding interface. 


Crystallization of SARS-CoV PLpro bound to the C-terminal domains of human and 


mouse ISG15. 


While X-ray crystal structures have been solved of SARS and MERS PLpro bound to Ub, 


our understanding of the interactions of these enzymes with ISGI5 has been largely limited to 


10 


enzymatic, mutational and computational modeling studies |" '° 


. To gain molecular insight into 
the specific interactions between PLPs and ISG15, we crystallized and determined the X-ray 
structures of SARS PLpro bound to the C-terminal domains of ISG15 from both human 
(ChISG15) and mouse (CmISG15). Attempts at co-crystallization of full-length ISG15s from 
human and other species with SARS CoV PLpro, as either various covalent adducts or non- 
covalently bound species, were made but were unsuccessful. However, the complex of SARS- 
CoV PLpro bound to the C-terminal domain of human ISG15 (ChISG15), which is the principle 
binding domain of ISG15, crystallized readily'’. Specifically, we utilized a form of ChISG15 
and CmISG15 that were modified with propargylamine at the c-terminus (ChISGI15-PA or 
CmISG15-PA) to form a suicide substrate that covalently links to the active site cysteine. From 


here on in, these covalently modified complexes of SARS PLpro are designated to as SARS- 


CoV PLpro-ChISG15 and SARS-CoV PLpro-CmISG15 for simplicity. 
X-ray SARS-CoV PLpro bound to the C-terminal domain of human ISG15. 


The X-ray structure of the SARS-CoV PLpro-ChISG15 complex was determined to 2.62 
A with two complete copies of SARS-CoV PLpro-ChISG15 within the asymmetric unit (Table 
3). The structure of CoV PLPs contains the classic tertiary fold associated with PLPs consisting 
of the finger, palm, thumb, and ubiquitin-like (Ubl) domains (Figure 4a). However, despite 
sharing the same tertiary fold, the manner in which SARS-CoV PLpro engages the C-terminal 
domain of hISGI5 is distinct from that of Ub. Overlaying the SARS-CoV PLpro-ChISG15 
structure with the structure of SARS-CoV PLpro with mono-Ub (PDB entry 4MM3) reveals that, 


compared to Ub, ChISG1S5 is shifted by approximately 12° in its global orientation (Figure 4b). 
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Closer examination of the binding interface uncovered key differences in how each 
substrate engages the surface of the SARS-CoV PLpro. Although there is only a minor variation 
surrounding a 180° flip of Trp123 between the two copies of ChISGI5 in the asymmetric unit, 
both ChISG15 chains are analogous in how they are accommodated by SARS-CoV PLpro 
compared to Ub (Figure 4e). The differences appear to be driven by how hISG15 and Ub differ 
in their ability to interact with a hydrophobic region consisting of residues Arg167, Met209, and 
Pro248. For Ub, the interactions with this region occur by means of a hydrophobic patch 
consisting of Ile44, Val70, and Leu8 'S In contrast, ISG15 lacks such a hydrophobic patch. The 
interaction is primarily driven by two distinct sets of hydrophobic interactions: hISG15’s Trp123 
with Met209 of the protease’s palm region, and hISG15’s Phel49 with a second hydrophobic 
site driven by Pro224 of the protease’s fingers region (Figures 4c, d). This predominate 
hydrophobic binding nature was also seen in the ITC results of SARS-CoV with ChSIG15, with 
an observed increase in entropic favorable interactions when compared with the hISG15 (Table 
2). Surprisingly, while the Ub-bound structure contains other polar and water-mediated 
interactions between Ub and the palm region of the PLpro, such interactions seem to be less 


pronounced in the more charged ChISGI5 structure Bs: 


Overall, there are only a few water 
molecules present within the interface between PLpro and ChISG15 compared to the larger 


number that are observed in the Ub-bound structure potentially indicating a lesser role of water 


mediated hydrogen bond interactions in the binding of hISGI5. 


To gain a better understanding of the residues mediating specificity for Ub or ISGI5, 
mutants targeting the palm and fingers regions of the of SARS-CoV PLpro were constructed and 
their catalytic activities towards various substrates were determined (Figure 5, Table 1). These 


substrates included Ub- and ISGI5-AMC, which are composed of the respective substrate 
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derivatized with the fluorescent leaving group 7-amido-4-methylcoumarin. Additionally, the 
small peptide Z-RLRGG representing the consensus recognition sequence for both Ub and 
ISG15 also attached to AMC was included to probe for potential changes affecting protein- 


protein interactions versus changes that impact intrinsic catalytic activity. 


One of the mutants, M209A, reduces the activity towards Ub-AMC by almost two-fold 
increases its activity towards Z-RLRGG-AMC by nearly five-fold, while retaining wildtype 
levels of activity for ISGI15-AMC. From a structural prospective, the reduction of Ub-AMC 
activity is not surprising as M209A shrinks the hydrophobic patch that engages the Ile44 
centered hydrophobic patch in Ub. The observed increase of the activity towards Z-RLRGG- 
AMC is less straightforward. M209 lacks direct interaction with the last five amino acids of Ub, 
or ISG15. However, its replacement by alanine could sterically open up a space that is more 
accommodating for the artificial Z-adduct of the peptide. The other two mutants were found to 
have increased catalytic efficiencies for processing either Ub or ISG15 with a corresponding 
reduction in activity for the other. R167E is over eight times more efficient than the wild type 
enzyme at hydrolyzing Ub but is about 20-times less efficient at hydrolyzing ISG15 (Figure 5c). 
In contrast, Q233E is nearly three-fold more efficient than wildtype at hydrolyzing ISG15 and 
two-fold less efficient hydrolyzing Ub. In the case of Ub, the mono-Ub structure suggests the 
charge flip with R167E may introduce an additional electrostatic interaction with either Arg42 or 
Gln49 of Ub. For the ChISG15 structure, this change occurs in close proximity to the interaction 
between Trp123 and Met209. Examination of the structure reveals Arg167 may contribute to z- 
m interactions involving Trp123 and Arg153 of ChISG15 (Figure 4e, Figure 5a). Replacing the 
longer and positively charged arginine with the shorter glutamate removes this contribution. 


Additionally, this charge flip may disrupt the electrostatic interaction between Arg153 and 
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Glul68 of the PLpro, potentially resulting in a loss of affinity for ChISG15. For Q233E, the 
change from a neutral polar to a charged group may create a potential electrostatic repulsion with 
the backbone carbonyl of Ala46 in Ub (Figure 5b). For the ChISG15 structure, in comparison, 
there is not a clear direct interaction between PLpro and ISG15 that would be affected. This 
suggests that the impact of this mutation may stem from internal changes within the PLpro 


regarding the flexibility of the finger region rather than direct interactions with the substrate. 
X-ray Structure of SARS-CoV PLpro bound to the C-terminal domain of mouse ISG15. 


In light of the differences in the nature of interactions between different species’ ISG15s 
revealed by ITC, as well as the apparent plasticity that can occur between Ub and hISG15 within 
the SARS-CoV PLpro active site, insights into the structural sources of this phenomenon were 
sought. To this end, a structure of SARS-CoV PLpro with the C-terminal domain of mISG1I5 
(CmISG15) was determined to a resolution of 2.4 A (Figure 6a). The structure of the catalytic 
core domain of SARS-CoV PLpro is consistent with the structure of the ChISGIS SARS-CoV 
PLpro structure. However, a major difference in the orientation of the Ubl domain of SARS-CoV 
PLpro is observed when bound to CmISG15. Unlike previous SARS-CoV PLpro structures 
where the Ubl domains differs only slightly in its position/orientation to the catalytic domain or 
is unobservable because of weak electron density suggesting potential flexibility, the Ubl in the 
SARS-CoV PLpro-CmISG15 structure takes a sharp, almost 90° turn in orientation when 
compared with previous X-ray structures (Figure 6c). The new orientation does not appear to 
form contacts with the bound CmISG15 suggesting that there is no direct influence on the Ubl to 
adopt this conformation. Further analysis of the new orientation reveals a seam made up of 
several electrostatic interactions between the Ubl domain and the thumb domains of the PLpro, 


indicating that this orientation can be stabilized (Figure 6b). 
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The Ubl domain was not the only global difference between the SARS-CoV PLpro- 
CmISG15 and SARS-CoV PLpro-ChISGI5 complexes. Surprisingly, the global orientation of 
CmISG15 is tilted 27° away from the fingers in SARS-CoV PLpro (Figure 6d). Interestingly, 
four amino acids appear to explain the lack of CmISG15’s accommodation in the same 
orientation to that of bound ChISGI5. At position 149 in ChISG15 and 147 in CmISGI5, there 
is an amino acid difference of a phenylalanine and an isoleucine respectively (Figure 6e). This 
change in CmISGI5 eliminates the hydrophobic interaction with the fingers region of SARS- 
CoV PLpro seen in the ChISG15 bound structure and likely aids in the tilt toward the Ubl 
domain. In addition, there is a change from Asn89, Thr125, and Asn151 in ChISG15 to Glu87, 
Ser123, and His149 in CmISG15 (Figure 6f). These amino acid differences create a hydrogen 
bond network between Glu87, His149, and Ser123 in the CmISGI15, which is absent in the 
ChISGI5 bound structure. This network locks the Glu87 into an unfavorable position for 
CmISG15 to bind in the same manner as ChISG15. The unfavorable position ultimately results 


in repulsion of CmIG15 against PLpro. 


In concert with a loss of the hydrophobic interaction with the finger region and potential 
electrostatic repulsions CmISG15 would incur in binding in the same orientation of ChISG15, 
the CmISG15 bound orientation facilitates the formation of numerous favorable interactions. 
Consistent with the X-ray structures of other de[SGylases bound to the C-terminal domain of 
human ISG15 7” ae the conserved Trp121 of mISGI15 is centric to the interface. However, 
unlike the interaction observed in the SARS-CoV PLpro-ChISGI5 structure, Trp121 in the 
mISG15 does not insert into a hydrophobic pocket. Instead, mISG15’s Trp121 and Pro128 form 
significantly smaller hydrophobic interaction with SARS-CoV PLpro via the proteses’s extended 


alkyl chain of Glul68 generated by the hydrogen bond formed between it and Thr171. In 
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contrast to the bound ChISG15 and outside the last five C-terminus consensus amino acids of Ub 
and ISGI5, this weak hydrophobic interaction between CmISGI5 and the protease is the only 
one observed. In addition to the hydrophobic binding contributions of mISG15 Trp121, a water 
molecule is observed to mediate a hydrogen-bonding network between Trp121’s indole nitrogen 
and several nearby PLpro residues (Figure 6g/h). Including this hydrogen bond network, 
CmISG15 forms almost 40% more hydrogen bonds with the protease then observed with 
ChISGI5. One set of these additional electrostatic interactions are centered around SARS-CoV 
PLpro’s Arg83 and CmISG15’s Lys132. Others are also observed to form between SARS-CoV 
PLpro’s Asn157, Gln175 and Argl67 and CmISG15’s Glu130, the carbonyl of Gly126, and 


His149 respectively (Figure 6g/h). 


Taken overall, the interface of bound CmISG15 of 872.5 A’ is similar in size to that of 
the ChISG15 ~767 A? interface, but consists of predominately electrostatic interactions whereas 
the former has a significantly greater hydrophobic component. In line with this, a greater 
number of water molecules, approximately 10, can be observed within the interface further 
indicating the degree to which binding may be driven more hydrophilic interactions. This 
characterization of the binding interaction is consistent with the thermodynamic properties 
observed for the binding of CmISG15 to SARS-CoV PLpro (Table 2). Also, when viewed in 
conjunction with the affinity data of mouse and human ISG15s as well as their C-terminal 
domains counterparts suggests that the type of interface formed between CmIS15 and SARS Co- 
V PLpro is less stable then its ChISG15 counterpart and may be one of the contributing factors to 
the overall reduction in binding affinity observed for full-length mISG15 when compared to 


hISG15 via ITC. 
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X-ray structure of the full-length, unbound form of mouse ISG15 and comparison human 


ISG15 


Whereas some viral deISGylating enzymes, such as vOTUs, are thought to exclusively 
interact with the C-terminus of ISG15, SARS PLpro has been shown to contain two Ub-binding 
sites that likely engage the N-terminus of ISG15s |" 'S “With only the full-length hISG15 
structure available, the impact of ISGI5 biodiversity within the N-terminal region has been 
difficult to assess. To address this issue, the X-ray structure of mISGI5 with both domains was 
determined to assess conformational variability of mISG15 compared to hISGI5. Initially, only 
a low resolution (4 A) structure could be determined. Truncation of the last five non-structured 
amino acids led to a structure that could be determined to a higher resolution of 3.25 A. Using 
the program Define Secondary Structure of Proteins (DSSP), the secondary structure of mISG15 
was determined to be comparable to hISG15. mISG15 contains nine beta sheets, two helices, but 
only two 39 helices within the C-terminal domain; this differs from the hISG15 which contains 
two 30 helices per domain of ISGI5 (Figure Fa): The asymmetric unit consists of 2 twisting 
filaments, containing 10 copies of mISG15. Upon examining the differences between the 
monomers within the asymmetric unit, some flexibility between N- and C-terminal domains was 
observed (Figure 7b/c). This flexibility was probed by aligning the C-terminal domains of each 
mISG15 monomer and measuring the angle between the point of divergence and the farthest N- 
termini. There is a 20.2° range of motion between the N- and C-terminal domains of mISG15 
(Figure 7b). Comparable flexibility was also observed in hISGI5 structures, with a 18.3° range 


of motion (Figure 7c) (PDB: 1Z2M, 3R66, 3PSE) *”*?*°. 


While each individual domain of hISG15 and mISG15 are similar in secondary structure 


and both show some structural variability between monomers, there is a drastic difference in the 
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overall conformations of the tertiary structures between mISGI5 and that of hISG15. This 
surprising difference is observed in the overall arrangement of the N- and C-terminal domains. 
Specifically, the twist about the C- and N-terminal domain of mISG15 in relation to hISG15 
ranges from 43.0° to 66.9° (Figure 8a). Closer investigation reveals that differences in tertiary 
arrangement can be attributed to several molecular interactions within the structures driven by 
the primary sequence differences of mISG15 compared to hISGI5. Specifically, the presence of 
Asp79 creates a kink in the hinge region of hISG15 as.a result of the carboxylate group of Asp79 
forming a hydrogen bond with the hydroxyl group of Thr101 (Figure 8b). This interaction does 


not occur in mISGI5 since Ser77 replaces Asp79. 


While amino acids in the hinge may play a part in the different conformational trajectory 
allowed in the region between the two domains in mISG15, the twisted feature of mISG15 may 
be stabilized by amino acid interactions in the core region of the protein. This core region is 
centered near Phe41, at the interface between the two domains. The hydrophobic interaction 
between the N-terminal and C-terminal domain at Leul34 and Phe41 may stabilize mISG15s’ 
twisted arrangement (Figure 8b). Other residues in close proximity, such as Pro39, appear to 
further stabilize this hydrophobic interaction in mISG15. Specifically, Pro39 in mISGI5 forms 
hydrophobic interactions with Phe41 contributing to the stabilization of mISGI15s’ tertiary 


arrangement (Figure 8b). 


In addition to the presence of favorable interactions, steric clashes are also likely 
responsible for preventing the occurrence of a shared tertiary arrangement between the mISG15 
and hISG15 species. When the N-terminal and C-terminal domains of mISG15 are configured to 
the arrangement of hISG15, there are no obvious clashes that prohibit mISG15 from adopting 


this conformation (Figure 8c). However, hISG1I5 may be incapable of adopting the mouse 
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conformation due to the potential steric clash between His39 and Glu139 that typically form a 
water mediated interaction (Figure 8c). Thus, the presence of the His39 residue may prevent 


hISG15 from configuring to the mISG15 conformation. 


Beyond the effects on tertiary structure, the sequence divergence between mISGI5 and 
hISG15 also impacts the potential binding surfaces. These differences in amino acid sequence 
also give rise to different electrostatic potential plots. Unlike mISG15, hISGI5 contains an 
uninterrupted negative surface band spanning across the two domains (Figure 8d). When the 
surface of mISG15 corresponding to the same area is examined, this negative band dissipates as 


the comparable area on mISGI5 is scattered with positive and neutral areas. 


DISCUSSION 


Differences in hISG15 and mISGI15 structure influence recognition by PLPs 


Among Ub and Ub-like proteins, ISG15 is unique for more than its divergent amino acid 
sequence. Similar to other Ub-like proteins, ISG15 has not been observed to form polymeric 
chains like Ub *!. However, ISG15 is the sole family member of ubiquitin-like modifiers that is 
comprised of two Ub-like folded domains. These two domains are tethered by a polypeptide 
linker that would suggest that the two domains can move freely and independent of each other. 
However, the four structures of hISG15 alone, or bound to viral proteins suggest the opposite 


(Figure 7c) Bh RS. 206 20 


Akin to polymeric Ub, ISG15 appears to have relatively limited 
conformations between the two domains. The structure of mISG15 (Figure 7b) furthers this 


assertion with one important caveat, one species ISG15s domain configuration may not be 


necessarily representative of others, or ISG15s in general. 
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In both hISG15 and mISG15 there is a hydrophobic interface mediated by a conserved 
phenylalanine (Phe41, see Figure 8). The highly-conserved nature of this phenylalanine and 
other surrounding residues suggests that ISG15s in general likely utilize hydrophobic forces to 
adhere the two domains together (Figure 1b, residues involved are boxed in blue). However, as 
shown with mISG15, the influence of the core on orientation of the domains can vary. In 
mISG15, residues forming the core solely dictate its inter-domain orientation, whereas in hISG15 
additional electrostatic interactions involving Glul39 and His39, and to a lesser extent Thr101 
and Asp79, appear to play an additional role to favor one orientation over another (Figure 8b,c). 
Intriguingly, the pairing of Glu139 and His39 is extremely unique to hISG15. Although Glu139 
or another acidic residue is well conserved at that position, His39 is typically a proline. In 
hISG15, it appears that this favorable electrostatic interaction may promote hISG15’s domain 
orientations. Multiple crystal structures reveal this conformation to be consistent despite 
differences in space groups and crystallization conditions, suggesting the observed tertiary 
structure a stable and likely preferred conformation (Figure 7c). In addition to favoring 
hISG15’s conformation, His39 may also act as a steric block preventing hISG15 from adopting a 


similar tertiary structure and conformation as that of mISG15. 


The observation of a potential steric block may result in MHV PLP2 not being able to 
recognize hISG15 which is why no cleavage of pro-hISG15 is observed (Figure 2). In contrast, 
no apparent steric hurdle is present for mISG15 in adopting a hISGI5 like inter-domain 
orientation which maybe why SARS and MERS PLpro are able to readily recognize mISG15. 
The lack of a steric hurdle might suggest that there are more allowable domain arrangements of 
mISGI5 than hISGI5 in solution that allow for a broader spectrum of PLPs to recognize 


mISG15. The potential uniqueness of hISG15’s structure may also fall in line with the inability 
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of other viral proteins, such as influenza NS1, to effectively engage ISG15s beyond those of 
human and primates °° Moreover, the distinctiveness of hISG15 may also fit into the recent 


assertions that hISG15 plays a divergent immune regulation role in humans than in other animals 
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Viral USP-Like PLPs accommodation of Ub versus ISG15s 


Ub interactions with proteins possessing an ubiquitin binding site have been observed to 
be reliant largely on the involvement of a hydrophobic patch comprised of several residues 
surrounding Ile44 'S. The absence of such a comparable patch in ISG1I5 as well as the more 
generally charged surface of hISGI5 naturally spurred speculation on mechanisms behind how 
PLPs may engage ISG15s. The SARS-CoV PLpro-ChISG15 complex reveals that in PLPs the 
binding interface may have evolved to recognize specific features of ISG15s from different 
species outside the five-amino acid, C-terminal sequence (Arg-Leu-Arg-Gly-Gly), they share 
with Ub. For hISG15, this includes not only the highly conserved Trp123 of ISG15s, but also 
specific interactions that are unique to hISG15. For instance, Arg87 in hISGI5 is often 
substituted for a lysine residue that is too short form any water mediated interactions with the 
SARS-CoV PLpro’s Asp230 and Ser222. Interestingly, Phel49, which forms a second 
hydrophobic interaction site with the protease’s Pro224, is one of the three amino acids recently 
implicated in species specificity among nairovirus VOTUs *°. Residue changes between hISG15 
and mISGI5 at the other two positions, Asn151 and Asn89, appear to impede the 
accommodation of mISGI5 in the active site of SARS-CoV PLpro compared to the more 
favorable hISGI5 orientation. This appears to advocate that differences at these three ISG15 
positions have a broader range of influence beyond only one class of viral proteases and could 


conceivable represent an evolutionary pressure that underlies part of ISG15 sequence diversity. 
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Also, SARS-CoV residues that interact with these residues as well as SARS-CoV PLpro 
deubiquitinating and deISGylating altering mutants, RI67E, M209A, and Q233E, may offer 
advantageous starting points for developing SARS-CoV PLpros with directed shifts in substrate 


specificities. 


Beyond the C-terminal domains interaction with PLPs, the N-terminus of human ISG15 
has been previous proposed to interact with a ridge helix spanning the conical PLP thumb 
domain with the Ubl domain to enhance affinity. Currently, no X-ray structures of an PLP with 
full-length ISG15 exist and previous computational models utilizing existing SARS-CoV PLpro 
structures bound with mono-Ub have had difficulty reconciling the significant gap that occurs 
between hISG15’s N-terminal domain and the protease when anchoring C-terminal hISG15 
domain on bound Ub. The SARS-CoV PLpro-CxISGI5 structures offer two synergistic 
explanations. First, the 6.3 A shift of ChISG15 relative to bound Ub translates the N-terminal 
domain towards the ridge helix (Figure 9a). In addition, the alternate Ubl conformation found in 
the SARS-CoV PLpro-CmISGI15 structure reveals that such a conformation translates the ridge 
helix up to 14 A toward the location of the ridge helix (Figure 9b/c). Interestingly, comparing 
this model to the X-ray structure of SARS-CoV PLpro bound to K48-linked di-Ub highlights that 
different facets of the Ub-like fold found in the N-terminal domain ISGI5 are likely involved in 
the interaction (Figure 9d). Specifically, in this model a triple serine repeat, Ser20-22, and 
Glu27 located in hISG15’s B2a3 loop and a3 respectively point toward the protease’s ridge 
helix creating a surface that is available to engaged by the bevy of charged and polar residues 
located on SAR-CoV PLpro’s ridge helix that has been previously implicated in binding’ : 
These potential electrostatic interactions may contribute to the four times larger enthalpy binding 


component observed in the interaction of full-length over the C-terminal domain of hISG15 


22 


when being accommodated by SARS-CoV PLpro. The impact of mISG15’s divergent domain 
orientation from hISG15 is also apparent. Initial molecular modeling of the hISG15 structure 
onto the CmISG15-bound SARS-CoV PLpro structure with CmISGI5 as an anchor reveals a 
steric clash with the ridge helix. This is the same for hISG15 drawn from any of its known X-ray 
structures. However, the divergent inter-domain orientation found in full-length mISGI5 
structure determined here allows for the mISG15’s N-terminal domain to fit unobstructed (Figure 
9e). The significant domain-domain orientations differences between hISG15 and mISG15 
results in a different facet of mISG15 facing the protease. Specifically, mISG15’s a1B3 loop 
comprises the surface oriented towards the ridge helix presenting a polar interface that ITC 
suggests may be involved in forming additional electrostatic interactions beyond those of the C- 
terminal domain. However, when comparing the ITC data for full-length and its CmISG15 
counterpart, these additional interactions might potentially come at an ethalphic cost. This may 
suggest that different ISG15 N-terminal domains may engage the protease’s ridge helix to 
differing degrees and thermodynamic characteristics. The need for PLPs not only to 
accommodate surface differences between species ISG15s, but also divergent inter-domain 
orientations highlight the difficulties for a CoV PLP to be active to all species ISG15s. This 
could also be perceived as a possible benefit to ISG15’s unique tandem Ub-like arrangement. 
Also, the improved structural perspective on how the N-terminus of ISG15s from different 
species may engage CoV PLPs provides a clearer path towards utilizing this region to influence 
deIGylating activities of CoV PLPs. As a result, combining this information along with the 
alterations possible in the C-terminal domain, fresh tools to addressing the role of 


deubiqutination and deISGylation through the use of reverse genetics systems can be envisioned. 


Possible evolution of CoV PLP recognition of species variances in ISG15s 
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Overall, the X-ray structural, enzymatic and biophysical data point to CoV PLP 
deISGylase activities being sensitive to species specific amino acid differences within ISGI5. 
Intriguingly, SARS-CoV and MERS-CoV PLpros, whose viruses replicate in a wide range of 
hosts, recognize and cleave proISG15 from almost all of the species tested. In contrast, the 
mouse specific MHV PLP2 is limited predominantly to the mouse substrate. This wider range of 
specificities appears substantially larger than that found recently in nairovirus VOTUs. This may 
imply that with the greater binding interface provided by the palm, fingers and thumb domains of 
the ubiquitin specific protease fold, CoV PLPs can either engage a wider array of ISG15s then 


that of vVOTUs, or be highly selective for just one or two ISG15s as in the case of MHV. 


Intriguingly, in all of the PLPs examined, their potential ability to engage different 
species ISGI5s is not fully restricted to only those from species their parent viruses infect. 
Understandably, this could be a by-product in the evolutionary process of a specific viral 
protease seeking to optimize towards a certain species ISGI5 and inadvertently picking up 
enzymatic activity towards another ISG15. Or, in an environment where a virus only has to 
optimize its replication in one species, there is less selective pressure and hence a greater chance 
of losing the ability of recognizing ISG15s from other species. Alternatively, some of these types 
of off-species PLP activities could be indicative of evolutionary memory for ISG15s. In essence, 
this could give a possible view into the zoonotic evolutionary history, or potential future 
zoonotic drift for a certain virus. Naturally, a wider sampling of CoV PLPs affinities for certain 
species ISG15s, knowledge of what species their parent viruses infect, and a multitude of reverse 
genetics experiments will be necessary to discern which of the above scenarios in taking place in 


the evolution of virus recognition of host ISGI5s. 


MATERIALS AND METHODS 
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Construct, Expression and Purification of PLPs for the ISG15 Protease Activity Assay and 
ITC. 


MHV PLP2 was expressed and purified as previously described “#33 SARS-CoV PLpro 
in expression vector pET2la and MERS-CoV PLpro in pET15b were transformed into 
Escherichia coli BL21(DE3) competent cells (New England Biolabs) by heat shock. Cells were 
grown at 37°C in LB broth supplemented with 100 ug/mL of ampicillin to OD¢o0 0.6-0.8 and 
expression induced with 0.8 mM isopropyl-f-D-thiogalactopyranoside (IPTG) at 25°C overnight 
for SARS-CoV PLpro and | mM IPTG at 18°C overnight for MERS-CoV PLpro. Cells were 
collected by centrifugation at 6000xg for 10 minutes and stored at -80°C. Cells were lysed in 
Buffer A (20 mM Tris [pH 7.5 for SARS-CoV PLpro, pH 7.0 for MERS-CoV PLpro], 500 mM 
NaCl, 10 mM f£-mercaptoethanol [BME]) supplemented with lysozyme for 30 minutes at 4°C, 
followed by sonication on ice at 50% power with a 50% duty cycle for a total of 6 minutes. 
Insoluble protein was removed by centrifugation at 70,600xg for 30 minutes and the supernatant 
filtered through a 0.80 um filter. The clarified supernatant was flowed over high density nickel 
agarose beads (Gold Biotechnology, Olivette, MO) pre-equilibrated with Buffer A. The column 
was washed with 10 column volumes of Buffer A supplemented with 30 mM imidazole, and the 
protein eluted with 10 column volumes of Buffer A supplemented with 300 mM imidazole. The 
PLpro was further purified by size exclusion chromatography using a Superdex 200 column (GE 
Healthcare, Pittsburgh, PA) pre-equilibrated with 100 mM NaCl, 5 mM HEPES [pH 7.5 for 


SARS-CoV PLpro, pH 7.0 for MERS-CoV PLpro] and 2 mM DTT. 
Purification of SARS-CoV PLpro for Complexation and Crystallization. 


Purification of SARS-CoV PLpro in expression vector pET11a for complexing with the 


C-terminal domain of human ISG15 (ChISG15) and mouse ISG15 (CmISG15) was adapted from 
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the previously described method ’. The cells were chemically lysed by resuspending them in 
150 mL of Buffer B (20 mM Tris [pH 7.5], 10 mM BME) and lysozyme and incubated at 4°C for 
30 minutes. The suspension was sonicated on ice at 50% power with a 50% duty cycle for a total 
of 6 minutes and centrifuged for 30 minutes at 40,900xg. The cell lysate was filtered and 
subjected to a 40% ammonium sulfate fractionation then centrifuged again for 30 minutes at 
40,900xg. The resulting pellet was resuspended in 250 mL of 1 M ammonium sulfate, 20 mM 
Tris [pH 7.5], and 10 mM BME and incubated at room temperature for 1 hour. The suspension 
was filtered and loaded onto a 50 mL Phenyl-Sepharose CL-4B column (GE Healthcare, 
Pittsburgh, PA) equilibrated with 1.5 M ammonium sulfate, 20 mM Tris [pH 7.5], and 10 mM 
BME. The protein was eluted using a 10-column-volume gradient to 100% Buffer B and washed 
with 2 additional column-volumes of 100% Buffer B. The fractions were pooled together and 
then diluted fivefold with Buffer B.. The protein was loaded onto a MonoQ 10/100 column (GE 
Healthcare, Pittsburgh, PA) equilibrated with Buffer B. The protein was eluted using a 10- 
column-volume gradient to 100% of a buffer composed of 0.5M NaCl, 20 mM Tris [pH 7.5], and 
10 mM BME with the initial flow-through off the column collected. The flow through was 
concentrated and put into dialysis in a 50 mM NaCl, 20 mM Tris [pH 7.5], and 10 mM BME 


buffer at 4°C overnight. 
Construct, Expression and Purification of proISG15s and mature ISG15s. 


ISG15s from human (Homo sapiens; Accession: AAH09507.1), mouse (Mus musculus; 
Accession: AABO2697.1), northern tree shrew (Tupaia belangeri; Accession: AFH66859.1), 
sheep (Ovis aries; Accession: AF152103.1), dromedary camel (Camelus dromedarius; 


Accession: XP_010997700.1), vesper bat (Myotis davidii; Accession: ELK23605.1), and 
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jackknife fish (Oplegnathus fasciatus; Accession: BAJ16365.1) in both pro and mature forms 


were prepared as described elsewhere ae 
ISG15 Protease Activity Assay. 


Activity assays of SARS-CoV PLpro, MERS-CoV PLpro, and MHV PLP2 with purified 
northern tree shrew proISGI5 (pro-nsISG15), sheep proISG15 (pro-shISG15), fish proISG15 
(pro-fISG15), mouse proISG15 (pro-mISG15), camel proISG15 (pro-cISG15), bat proISG15 
(pro-bISG15), and human proISGI5 (pro-hISG15) were adapted from the previously described 


methods *. 
ITC of ISG15 with PLpros from MERS-CoV and SARS-CoV. 


ITC was performed using a Microcal PEAQ-ITC (Malvern, Worcestershire, UK). There 
were 19 injections of 2 uL each at 25°C with a reference power of 6 pcal/s. The mature forms of 
ISG15s along with PLpros from MERS-CoV and SARS-CoV were dialyzed at 4°C in 50 mM 
HEPES [pH 7.4], 200 mM NaCl, and 1 mM DTT. All experiments were run in duplicate. For 
direct binding experiments, 227 uM and 276 uM of SARS-CoV and MERS-CoV PLpro 
respectively were placed in the cell with 2.3-2.6 mM of mature hISGI5 in the syringe. For direct 
binding experiments of SARS-CoV with ChISG15, 303 uM was placed in the cell with 3.4 mM 
in the syringe respectively and for SARS-CoV with CmISG15, 393 uM was placed in the cell 
with 3.8 mM in the syringe respectively. For competitive experiments related to SARS-CoV 
PLpro, mixtures containing 100 uM of protease with 50 uM of either sheep ISG15 (shISG15), or 
mISGI5 was placed in the cell with 1 mM of mature hISGI5 in the syringe. For competitive 
experiments related to MERS-CoV PLpro, mixtures containing 220 uM and 270 uM of protease 


with 110 uM and 170 uM of shISGI5 and mISG15 respectively was placed in the cell. The 
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syringe contained 2.3 mM and 2.7 mM of mature hISGI5 in the syringe for shISG15 and 
mISGI5 assays respectively. The data were processed using Microcal PEAQ-ITC Analysis 


Software. 
Functional studies of SARS-CoV PLpro Mutants. 


The SARS-CoV pET-15b-PLpro mutants (residues 1541-1855 of the SARS-CoV viral 
polyprotein) were generated using site-directed mutagenesis and the QuickChange® approach 
(Agilent). Expression and purification for the wild type and each mutant of the SARS-CoV 


PLpro were performed as previously described ™, 


The steady-state kinetic parameters of SARS-CoV PLpro wild-type and mutants were 
determined for three different ubiquitin-based fluorescent substrates, utilizing 7-amino-4- 
methylcoumarin (AMC), commonly used to assess the protease, deubiquitinating, and 
deISGylating activity of PLPs, including a small peptide substrate, Z-RLRGG-AMC (Bachem), 
Ub-AMC (LifeSensors, Inc.), and ISGI15-AMC (Boston Biochem/R&D Systems). Kinetic 
assays with Ub-AMC and ISG15-AMC were performed on the same day and side-by-side in the 
same assay plate to directly compare the enzymatic activity of SARS-CoV PLpro to that of each 
of the mutants. The steady-state kinetic studies were also repeated for the wild-type and mutants 
approximately five months apart and the resulting duplicate data were combined for analysis. 
Kinetic assays with the peptide substrate were also performed in triplicate. For all experiments, 
the assay conditions, i.e. buffering conditions and assay volume etc, were setup as previously 
described **. The exception was that the stock substrates purchased from the vendors had 


different lot numbers. The steady-state kinetic data obtained from separate experiments 
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performed on different days and with different substrates lot numbers helped to ensure that the 


trends in the resulting kinetic parameters were reproducible. 


The enzymatic activity of PLpro-mediated hydrolysis of the fluorophore, AMC group 
was determined using a BioTEK Synergy H1 multimode microplate reader at 25°C with a 
wavelength of excitation at 360 nm (bandwidth=40 nm) and an emission wavelength of 460 nm 
(bandwidth=40 nm). The change in the relative fluorescence as a function of time (RFU/min) 
was monitored over a sufficient time period to allow determination of the enzymatic rate in the 
steady-state region. For the ISG15-AMC assay, the substrate concentrations were varied from 
0.2 uM up to 19.2 uM. The reactions were initiated by the addition of enzyme with the final 
enzyme concentrations as follows: 0.48 nM WT, 0.23 nM Q233E, 0.23 nM M209A, or 7.3 nM 
R167E. For the Ub-AMC assay, substrate concentrations were varied from 0.5 uM to 17.6 uM. 
The final enzyme concentrations were 3.7 nM WT, 7.3 nM Q233E, 7.3 nM M209A, or 0.23 nM 
R167E. For the Z-RLRGG-AMC assay, the concentrations of substrate were varied from 0.8 uM 
to 50 uM and the final concentration of the wild type enzyme was 0.14 uM. To capture the 
initial rate of peptide hydrolysis for the M209A mutant a lower enzyme concentration of 25 nM 
was used. As is typical for SARS-CoV PLpro, the enzyme could not be saturated with the Ub- 
AMC and Z-RLRGG-AMC substrates. As such, the kinetic response of the enzyme to these 
substrates was linear and thus the data were fit to a line to approximate the catalytic efficiency 
(Kea/Km) for each enzyme. For the ISG15-AMC assays, the data were fit to the Michaelis- 
Menten equation to determine the associated kinetic parameters (Kea, Km, and Kea/Km) for each 
enzyme ** Saturation was not attained with the R167E mutant enzyme for ISG15-AMC and 
therefore this kinetic data was also fit to a line to determine the apparent (Kea/K). The errors 


associated with each kinetic parameter were obtained from the best-fit line or curves for each 
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mutant. All data, from separate experiments, were included in the fits to arrive at the final errors 


(Table 1). 
SARS-CoV PLpro-CmISG15 and SARS-CoV PLpro-ChISG15 Complex Formation. 


Expression of CmISG15, or ChISG1I5 occurred using a vector pT YB2 and was purified 
as previously described to form a propargylamine-derivatized thioester product (CmISG15-PA, 
ChISG15-PA) ”°. Briefly, to obtain complex, purified protease was added directly to the mixture 
in equi-molar ratios, and incubated for 2-4 hours at room temperature and left 4°C overnight. To 
further purify the complex, anion exchange chromatography was used, eluting from a MonoQ 
10/100 column using a linear gradient from 0 to 1 M NaCl (SARS-CoV PLpro-CmISG15) or 
250 mM NaCl (SARS-CoV PLpro-ChISG15) with 50 mM Tris (pH 8.0 for SARS-CoV PLpro- 
CmISG1I5 and pH 9.0 for SARS-CoV PLpro-ChISGI5), followed by size exclusion 
chromatography on a Superdex 75 column (GE Healthcare, Pittsburgh, PA) pre-equilibrated with 
100 mM NaCl, 5 mM HEPES [pH 7.5], and 2 mM DTT. For SARS-CoV PLpro-ChISG15, an 
additional purification step prior to anion exchange chromatography was used to eliminate 
residual ChISGI5 by size exclusion chromatography using a Superdex 200 column (GE 
Healthcare, Pittsburgh, PA) pre-equilibrated with 100 mM NaCl, 50 mM Tris [pH 7.5], 2 mM 


DTT. 


Crystallization of SARS-CoV PLpro-CmISG15, SARS-CoV PLpro-ChISG15, and 


mISG15. 


Purified SARS-CoV PLpro-CmISG15, SARS-CoV PLpro-ChISG15, and mISG15 were 
screened against a series of Qiagen NeXtal suites by hanging drop using a TTP Labtech 


Mosquito (TTP Labtech, Herfordshire, UK) at 8.8 mg/ml, 8.88 mg/ml, and 16 mg/ml 
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respectively. For SARS-CoV PLpro-CmISG15, the initial screen yielded the best crystals in a 
solution containing 65% (vol/vol) MPD and 0.1 M Tris (pH 8.0). These crystals were optimized 
using the Additive HT Screen from Hampton Research. The final SARS-CoV PLpro-CmISG15 
crystals were obtained through vapor diffusion using a 500 wL reservoir with 4 uL hanging drops 
mixed 1:1 with protein solution and reservoir solution, which also contained 0.25 ul of 30% 
(w/v) Trimethylamine N-oxide dihydrate. For SARS-CoV PLpro-ChISG15, the initial screen 
yielded the best crystals in a solution containing 0.2 M lithium sulfate, 0.1 M Bis-Tris (pH 6.5) 
and 25% (w/v) PEG3350. The initial crystal conditions for SARS-CoV PLpro-ChISG15 crystals 
were optimized along salt and PEG3350 gradients. as well as using the Additive Screen from 
Hampton Research. The final SARS-CoV PLpro-ChISG15 crystals were obtained through vapor 
diffusion using a 500 uL reservoir with 4 wL hanging drops mixed 1:1 with protein solution (6.99 
mg/mL) and reservoir solution (0.1 M lithium sulfate, 0.1 M Bis-Tris, [pH 6.5], 22% PEG3350), 
which also contained 0.25 ul of 30% (v/v) glycerol. For mISG15, the initial screen yielded the 
best crystals in a solution containing 0.2 M ammonium sulfate, 0.1 M tri-sodium citrate pH 5.6, 
15% (w/v) PEG4000. This condition was further optimized along buffer, pH, and PEG4000 
gradients in addition to using the Additive Screen from Hampton Research. In conjunction with 
using these optimization methods, the mISG15 was shorted by five amino acids. The final 
crystals of the five amino acid shortened mISGI5 were obtained through vapor diffusion using a 
500 uL reservoir with 4 wL hanging drops mixed 1:1 with protein solution (6.9 mg/mL) and 
reservoir solution (0.2 M ammonium sulfate, 0.1 M sodium acetate pH 4.6, 12% PEG4000 and 


0.2 M sodium malonate). 


Crystals of SARS-CoV PLpro-CmISG15, SARS-CoV PLpro-ChISG15, and mISG15 


were collected and flash frozen in liquid Nz. Cryogenic solutions for SARS-CoV PLpro- 
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CmISG15, SARS-CoV PLpro-ChISG15 constituted their respective mother liquors. For 
mISGI5 crystals, they were passed from a 5% to a 12% solution 1:1:1 of glycerol, dimethyl 
sulphoxide (DMSO), and polyethylene glycol known as EDG >. Data sets were collected at the 
Advanced Photon Source (Argonne National Labs, Argonne, IL). A data set for SARS-CoV 
PLpro-CmISG15 was collected at the LS-CAT beamline 21G at a wavelength of 0.9786 A using 
a MAR300 detector, whereas data sets for SARS-CoV PLpro-ChISG15 and mISG15 were 
collected at the SER-CAT beamlines 22ID and 22BM at 1 A using MAR300hs detectors. All 


data sets were collected at 100 K. 
Data Processing and Structure Solutions. 


Data sets were indexed, integrated and scaled using HKL-2000 %6 All the structures were 
solved by molecular replacement using Phaser ie Subsequently, each structure was rebuilt 
initially using Autobuild *8 followed by successive rounds of manual model building and 
refinement using Coot °° and Phenix “°. The initial solution for the SARS-CoV PLpro-CmISG15 
complex was achieved by using the catalytic core of a previous SARS-CoV PLpro structure 
(PDB entry 3E9S). Density from the last 10 amino acids of the CmISG15 molecule served as an 
anchor for the initial global placement of the CmISGI5 from a Erve nairovirus VOTU-CmISG15 
complex (PDB entry 5JZE). This partial model was then used as a search model along with the 
Ub-like SARS-CoV PLpro domain from PDB entry 4MM3 to obtain a complete global model 
using Phaser °7 For the SARS-CoV PLpro-ChISG15 complex, an initial molecular replacement 
solution was obtained by using the core and other elements from the SARS-CoV PLpro- 
CmISG15 as a search model. A partial molecular replacement solution using Phaser *’ for 
mISGI5 was obtained by searching with the CmISG15 from the Erve nairovirus vOTU- 


CmISG15 complex structure (PDB entry 5JZE). This partial model was used in a sequential 
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Phaser *” run using a mISG1I5 N-terminal domain homology model based on the previously 
solved hISG15 structure (PDB entry 1Z2M). All structures were validated using Molprobity a 
and have good Ramachandran statistics: SARS-CoV PLpro-CmISG15 (96% favored and 4% 
allowed), SARS-CoV PLpro-ChISG1I5 (96% favored and 4% allowed) and mISG15 (99.08% 
favored and 0.92% allowed). All structures have been deposited in the protein data bank. Codes 


can be found on Table 3. 
Electropotential plots. 


Figure renderings involving electropotential plots were performed using the PDB2PQR 


server and the surface generated using the adaptive Poisson-Boltzmann solver (APBS) aesA 
ACCESSION NUMBERS 


All structures have been deposited in the Protein Data Bank (PDB code 5TL7 for SARS-CoV 
PLpro-CmISG15, 5TL6 for SARS-CoV PLpro-ChISG15, and STLA for mISG15). Source sequences for 
ISG15s were GenBank numbers AAH09507.1 for Homo sapiens, AAB02697.1 for Mus musculus, 
AFH66859.1 for Tupaia belangeri, ELK23605.1 for Myotis davidii, BAJ16365.1 for Oplegnathus 
fasciatus, and NCBI Reference Sequence XP_010997700.1 for Camelus dromedarius. Source sequences 
for PLPs were GenBank numbers AFS88944.1 for MERS CoV PLpro, AKP80587.1 for PEDV PLP2, 
AHM88399.1 for PDCoV PLpro, and UniProtKB/Swiss-Prot numbers PO0C6U8 for SARS CoV PLpro 


and POC6V0 for MHV PLP2. 
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Figure and Table Legends 


Table 1. Kinetic Parameters of SARS-CoV PLpro WT and mutants with different ubiquitin- 


based fluorescent substrates. 


Table 2. Isothermal Titration Calorimetry. 


Table 3. Data collection and refinement statistics. 

Figure 1 (a) The PLPs from SARS CoV (Accession number POC6U8), MERS CoV (Accession 
number AFS88944.1), MHV (Accession number POC6V0), PEDV (Accession number 
AKP80587.1), and Porcine deltacoronavirus (PDCoV) (Accession number AHM88399.1). The 
sequence numbering ruler is based on the SARS-CoV PLpro sequence. Residues forming the 
catalytic triad are boxed in purple. The secondary structure of SARS-CoV PLpro based on 
DSSP is denoted, with a- and 3;0-helices in reddish orange and B-sheets in wheat. Residues 
forming the “ridge” of SARS-CoV PLpro are marked with an asterisk (*). Residues of SARS- 
CoV PLpro mutated in this study are denoted (“x”). (b) Sequence alignment of ISG15s from 
mouse (Mus musculus; Accession: AABOQ2697.1), human (Homo sapiens; Accession: 
AAH09507.1), dromedary camel Camelus dromedarius; Accession: XP_010997700.1), vesper 
bat (Myotis davidii; Accession: ELK23605.1), sheep (Ovis aries; Accession: AF152103.1), 
northern tree shrew (Tupaia belangeri; Accession: AFH66859.1), and jackknife fish 
(Oplegnathus fasciatus; Accession: BAJ16365.1) along with ubiquitin. The sequence numbering 
ruler is based on mouse ISGI5. Secondary structure of mouse ISG15 based on DSSP is denoted 
in purple. The hinge region between B4 and £5 is denoted using a dashed line in teal. 
Percentages indicate sequence identity relative to mouse ISG15. Residue positions identified in 


this study that form key interactions with SARS-CoV PLpro are boxed in red. Residue positions 
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implicated in driving the specific tertiary fold features of mouse and human ISGI5 are boxed in 


blue. 


Figure 2 Cleavage profiles of SARS-CoV PLpro, MERS-CoV PLpro, and MHV PLP2 with 
proISG15s derived from human, sheep, northern tree shrew, jackknife fish, mouse, dromedary 
camel, and vesper bat. Reactions were conducted at 37°C with 10 uM of each proISG15 and 20 


nM of each PLpro/PLP2 with samples taken at the indicated time points. 


Figure 3 ITC binding isotherms with the raw heat (top panel) and the integrated heats of 


injection (bottom panel) shown for each interaction. 


Figure 4 (a) Cartoon representation of the SARS-CoV PLpro complexed with ChISG15 (blue). 
Secondary structure of the PLpro is designated, with helices and loops rendered in reddish 
orange and B-sheets rendered in. wheat. The structural domains of the PLpro are identified, 
consisting of the fingers (white), palm (purple), thumb (green), and Ubl (red) labeled with 
colored boxes. The vinyl thioether propargylamine linker is colored green, and the Zn(ID ion in 
light purple. (b) Overlay of ChISG15 (blue) bound to SARS PLpro (reddish orange) compared 
to a Ub (light brown)-bound structure (PDB entry 4MM3) based on a secondary structure 
alignment of the respective PLpros. The approximate degree of shift in orientation of analogous 
a-helices is indicated. (c) Comparison of the interaction of ChISGI5 (blue) versus Ub (light 
brown) with the hydrophobic patch of SARS-CoV PLpro, with the site of interaction within the 
overall structure indicated by an arrow. Side chain conformations of the PLpro are colored 
orange for the ChISG15 bound structure, gray for the Ub-bound structure. Residues E127-D133 
of ChISG1I5 and A46-T55 of Ub were removed for clarity. Inter- and intramolecular distances 


(A) are shown by black dashes. (d) Comparison of ChISG15 versus Ub at the site of the 
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additional hydrophobic interaction in the ChISG15-bound structure. Colored as in (c). 
Intermolecular distances (A) are shown as as black dashes. (e) Interchain variability of ChISG15 
Trp123 in binding to SARS-CoV PLpro. The two copies within the asymmetric unit were 
overlaid based on a secondary structure alignment of the PLpros. Side chains of the PLpro 
colored in orange correspond to the ChISG1I5 colored light blue, and the PLpro side chains 
colored burgundy to the ChISG15 colored a darker blue. Intermolecular distances are colored in 


black. Residues E127-D133 of ChISG15 omitted for clarity. 


Figure 5 (a) and (b) Sites within the ChISG15- and Ub-bound SARS-CoV PLpro structures 
corresponding to the mutations causing increased Ub-AMC activity (top) versus increased 
ISG15-AMC activity (bottom). Structures shown as cartoons, colored as in (4c). Intermolecular 
distances to indicate proximity for SARS-CoV PLpro-Ub are shown as black dashes with the 
ones for SARS-CoV PLpro-ChISGI5 in purple. (c) Activity of SARS-CoV PLpro mutants 


towards Ub-, and ISG15-AMC. Corresponding data found in Table 1. 


Figure 6 (a) Cartoon representation of the SARS-CoV PLpro complexed with CmISG15 (green). 
Secondary structure of the PLpro is designated, with helices and loops rendered in gold and B- 
sheets rendered in silver. The structural domains of the PLpro are identified, consisting of the 
fingers (white), palm (purple), thumb (green), and Ubl (red) labeled with colored boxes. (b) 
Cartoon representation in wall-eye stereo view, with transparent surfaces, of the interactions 
between the thumb domain (gold) and Ubl-domain (gray) for the SARS-CoV PLpro-CmISG15 
bound structure. Amino acids involved in interactions are represented as sticks and hydrogen 
bonds are indicated by black dashes and distances (A) labeled in red. (c) Ribbon overlay of 
SARS-CoV PLpro complexed with CmISG15 (gold) and SARS-CoV PLpro complexed with 


ChISGI5 (reddish orange) demonstrating the movement of the Ubl-domain. (d) Overlay of 
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SARS-CoV PLpro (gold) complexed with CmISGI5 (green) and the ChISG15 (blue) from the 
SARS-CoV PLpro complex based on a secondary structure alignment of the respective PLpros. 
The approximate degree of shift in orientation of analogous a-helices is indicated. (e/f) 
Comparison between the SARS-CoV PLpro (reddish orange) complexed with ChISGI5 (blue) 
and CmISG15 (green). Electrostatic interactions indicated with black dashes and distances (A) 
labeled in red or yellow, respectively. (g/h) Interactions between SARS-CoV PLpro (gold) and 
CmISGI15 (green) with electrostatic interactions indicated with black dashes and distances (A) 


labeled in red and hydrophobic interactions indicated by the yellow transparent surface. 


Figure 7 (a) Secondary structure of monomer A of full-length mISG15 (purple) with B-sheets 
colored in silver. (b) Monomers of mISG15 were overlaid, aligning the C-terminal domains, to 
observe the flexibility of the monomeric units. (c) Monomers of hISGI5 (green; PDB 1Z2M, 
gray; PDB 3R66, blue; PDB 3PSE) were overlaid, aligning the C-terminal domains, to observe 


the flexibility of the monomeric units. 


Figure 8 (a) Overlay of mISGI5 (purple) and hISG1I5 (green; PDB 1Z2M) shown in ribbons. (b) 
Comparison of interactions responsible for the structural conformations of hISG15 (green) and 
subunit A mISG15 (purple) with dashed lines representing bond lengths in A. (c) Interactions 
observed when mISG15 is forced into the hISG15 conformation (left panel) and vice versa (right 


panel). (d) Electrostatic surface of mISG15 and hISG15. 


Figure 9 (a) SARS-CoV PLpro (reddish orange) from the ChISG15 bound structure with full- 
length human (green; PDB 1Z2M) overlaid with full-length human based off of a secondary 
alignment using Ub (light gray; PDB 4MM3). The distance between the two-overlaid full-length 


human structures is measured in A and is indicated by a blue dashed line. For all figures the 
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highlighted yellow portion demonstrates the secondary binding site ridgeline. (b) SARS-CoV 
PLpro (gold) from the CmISG15 bound structure with full-length human (green; PDB 1Z2M) 
based off of a secondary structure alignment using ChISGI5. (c) Close-up overlay of the 
secondary binding site helix from SARS-CoV PLpro (gold) from the CmISG15 bound structure, 
SARS-CoV PLpro (reddish orange) from the ChISG15 bound structure, and SARS-CoV PLpro 
(dark gray; PDB 5E6J) from the K48 di-Ub bound structure. The distance between the helix 
from SARS-CoV PLpro (gold) from the CmISG15 bound structure and the helix from SARS- 
CoV PLpro (reddish orange) from the ChISG15 bound structure is measured in A and is 
indicated by a blue dashed line. (d) SARS-CoV PLpro (reddish orange) from the ChISG15 
bound structure with full-length human (green; PDB 1Z2M) overlayed with K48 di-Ub (dark 
brown; PDB 5E6J) both based off of a secondary alignment using ChISGI5. (e) SARS-CoV 
PLpro (gold) from the CmISGI5 bound structure with full-length mouse overlayed with full- 
length human (green; PDB 1Z2M) both based off a secondary structure alignment using 


CmISG15. 
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Figure 9 
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SARS-CoV PLpro 


Graphical abstract 


Table 1. Kinetic Parameters of SARS PLpro WT and mutants with different ubiquitin-based fluorescent 


substrates. 


Substrate SARS PLpro Enzymes 


Kinetic Parameter WT R167E M209A Q233E 


RLRGG-AMC’ 


0.062 + 0.001 0.727 + 0.007 0.211 + 0.002 


Kea/Kim (UM min’') 0.141 + 0.002 
Ub-AMC? 
Kea/Kim (UM min") 7.22 + 0.56 58.0 £2.5 2.90 + 0.32 3.14 40.32 
ISG15-AMC” 
Kea/Kim (UM! min”) 50.7 + 9.0 2.76 + 0.10" 54.3 + 12.8 132 +23 


“Best-fit slope values derived k,,,, for nonsaturating substrates and approximates kea/Km 


Steady-state values following Michaelis-Menten kinetics are reported from duplicate measurements 


54 


Table 2. Isothermal Titration Calorimetry. 


Protein 


SARS-CoV PLpro 


Human 
C-Human 
Sheep 
Mouse 


C-Mouse 


MERS-CoV PLpro 


Human 
Sheep 


Mouse 


N 


(sites) 


0.932 + 0.032 
1.31 40.015 
1 
1 


1 


0.785 + 0.027 
1 


1 


Kp 


(uM) 


20.50 + 4.48 
57.6 + 3.21 
200 + 41 
196 + 68 


1870 + 173 


59.3 +12.7 
172 + 38.5 


376 + 52.8 


AH 


(kJ/mol) 


-27.20 + 1.90 
-4.03 +.019 
-106 + 10 
-92 + 15 


-20.71 + 0.486 


-11.34 + 0.268 
-19.92 + 0.420 


-43.51 + 0.535 


AG 


(kJ/mol) 


-26.80 
-24.18 
-21.10 
-21.20 


-15.56 


-24.15 
-21.51 


-19.54 


-TAS 


(kJ/mol) 


Table 3. Data collection and refinement statistics. 


Data collection 
Space group 
Cell dimensions 

a, b, c (A) 

a, B,y ©) 
Resolution (A) 
Ryym OF Rmerge 
I/ol 
Completeness (%) 


Redundancy 


Refinement 
Resolution (A) 
No. reflections 
Rwork (%)/ Riree (%) 
No. atoms 
Protein 
Ligand/ion 
Water 
B-factors 
Protein 
Ligand/ion 
Water 


R.m.s. deviations 


SARS-CoV 


PLpro-ChISG15 (PDB 
entry STL6)* 


P2);2)2; 


46.9, 87.0, 221.5 
90, 90, 90 

50.00-2.62 (2.67-2.62)' 
11.1 (58.3) 

12.7 (2.07) 

98.4 (96.0) 


5.0 (3.9) 


38.67-2.62 (2.71-2.62) 
27,644 


19.1/25.0 


6177 
15 


133 


49.27 
59.94 


41.95 


SARS-CoV 


PLpro-CmISG15 (PDB 
entry 5TL7)* 


P2)2)2) 


76.0, 98.2, 106.6 

90, 90, 90 

50:00-2.44 (2.53-2.44)' 
9.9 (74.2) 

18.0 (2.71) 

99.1 (99.8) 


2.4 (2.1) 


32.20-2.44 (2.53-2.44) 
29,886 


18.7/25.3 


6218 
10 


231 


61.13 
57.02 


58.23 


56 


mISG15 (PDB entry 
5TLA)* 


P2 


85.8 60.2 172.7 
90, 93.1, 90 

50.00-3.25 (3.31-3.25)" 
5273) 

17.4 (2.27) 

98.4 (88.0) 


2.4 (2.1) 


42.18-3.25 (3.36-3.25) 
27,929 


25.6/30.9 
11,491 
0 


0 


86.7 


Bond lengths (A) 0.004 0.008 0.002 
Bond angles (°) 0.54 1.1 0.57 


*Datasets collected from a single crystal for each structure. ‘Values in parentheses are for highest-resolution shell. 
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Highlights 


Coronavirus PLPs possess immune modulating deubiquitinating/deISGylating activity 
Different species’ ISG15s have different binding characteristics with PLPs 

New structures of PLPs with an ISG15 domain expose different binding orientations 
A structure of full length mouse ISG15 reveals overall differences from human ISG15 
ISG15 variation affects PLP function and potentially host immune responses 
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